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Abstract 

Background: The genus Limonium Miller comprises annual and perennial halophytes that can produce sexual 
and/or asexual seeds (apomixis). Genetic and epigenetic (DNA methylation) variation patterns were investigated in 
populations of three phenotypically similar putative sexual diploid species (L nydeggeri, L ovalifolium, L lanceolatum), 
one sexual tetraploid species (L vulgore) and two apomict tetraploid species thought to be related (L dodartii, 
L multiflorum). The extent of morphological differentiation between these species was assessed using ten diagnostic 
morphometric characters. 

Results: A discriminant analysis using the morphometric variables reliably assigns individuals into their respective 
species groups. We found that only modest genetic and epigenetic differentiation was revealed between species by 
Methylation Sensitive Amplification Polymorphism (MSAP). However, whilst there was little separation possible 
between ploidy levels on the basis of genetic profiles, there was clear and pronounced interploidy discrimination on 
the basis of epigenetic profiles. Here we investigate the relative contribution of genetic and epigenetic factors in 
explaining the complex phenotypic variability seen in problematic taxonomic groups such as Limonium that operate 
both apomixis and sexual modes of reproduction. 

Conclusions: Our results suggest that epigenetic variation might be one of the drivers of the phenotypic divergence 
between diploid and tetraploid taxo and discuss that intergenome silencing offers a plausible mechanistic explanation 
for the observed phenotypic divergence between these microspecies. These results also suggest that epigenetic 
profiling offer an additional tool to infer ploidy level in stored specimens and that stable epigenetic change may play 
an important role in apomict evolution and species recognition. 
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Background 

DNA sequence divergence clearly plays a leading role in 
shaping the phenotypic variation observed between most 
taxa e.g. [1,2] but does not explain all forms of adaptive 
phenotypic differentiation [3-5]. Epigenetic modifications 
of DNA and histones, the core components of chromatin 
are able to influence the expression of the underlying 
genes and so phenotype [6]. Of these epigenetic mecha- 
nisms, the best studied and the one to show more ex- 
amples of transgenerational stability is cytosine DNA 
methylation [6,7]. Evidence has also been found sug- 
gesting that heritable phenotypic variation observed in 
natural populations can be due to stable epigenetic 
variation, even in the absence of genetic variation and 
could play a role in plant adaptation and evolution 
[5,8-10]. The scope for epigenetic mechanisms driving 
morphological differentiation is perhaps best illustrated 
when genetic variation is lacking. Understanding the 
relative importance of genetic and epigenetic sources of 
phenotypic variation where both systems are operating 
is therefore attracting increasing interest. For example, 
Laguncularia racemosa is a mangrove plant species that 
shows low genetic variability but in populations from 
distinct, nearby habitats, cytosine methylation variation 
among individuals correlates more closely with envir- 
onmental variation than does genetic variation [11]. In 
the perennial Viola cazorlensis, cytosine methylation pat- 
terns were found to be partitioned and positively correlated 
with adaptive genetic variation [12]. Also, in populations 
of individuals with reduced or negligible genetic vari- 
ation such as those of triploid asexual dandelion lineages 
(apomixis; diplospory), changes in genomic methylation 
patterns are found between individuals [13]. 

The genus Limonium Miller (sea-lavenders; Plumbagi- 
naceae) has long been recognized to have a history of 
recurrent hybridization and polyploidization, and com- 
prises 150 [14] to 350 taxa recognized across coastal, 
steppe and desert regions (e.g. [15,16]). This wide range 
is due to the description of new taxa, mainly micro- 
species from geographically restricted areas. In this genus, 
a sporophytic self-incompatibility system is linked with 
pollen-stigma dimorphisms, A-pollen type grains ger- 
minate on papillose stigmas and B-pollen type germinate 
in cob-like stigmas, while the complementary combina- 
tions produce no successful fertilization [17,18]. Most 
sexual species of Limonium usually have a dimorphic 
self-incompatibility system (both pollen and stigmas 
are dimorphic) while agamospermous species are gen- 
erally monomorphic and have monomorphic populations 
[14,19]. Determination of these characters in individuals 
from natural and/or experimental populations has since 
long been used as an indirect method for estimation of 
each species reproduction mode [14,17,18]. The high 
number of polyploid taxa has been explained to be a 



consequence of this self-incompatibility system and the 
ability of polyploid hybrids to produce seeds asexually 
via apomixis (agamospermy; asexual seed formation) 
[14,17-20]. Ixeris-type embryo sacs with non-haploid eggs 
are found in triploid (2n = 3x = 27) Statice oleaefolia var. 
confusa [20]. In triploid and tetraploid agamospermous 
species of the L. binervosum (G. E. Sm.) Salmon group 
diplospory followed by parthenogenesis is reported [21,22]. 
Molecular phylogenetic studies have tried to resolve the 
taxonomic complexity within this genus in a global per- 
spective using nuclear DNA sequence information [23] 
and plastid DNA [24-26]. 

In Continental Portugal about 15 Limonium species 
have been recognized with ecological importance for plant 
communities of the Atlantic and Mediterranean coastlines 
[15,27]. Among these, the L. ovalifolium complex consists 
a group of three sexual diploids (2n = 2x =16): L ovalifolium 
(Poir.) O. Kuntze, L. nydeggeri Erben and L. lanceolatum 
(Hoffmanns & Link) Franco [28]. The first species 
has a broader distribution including several Sites of 
Community Importance (SCI) for the Mediterranean 
biogeographical regions [29] in the West (Estremadura), 
Southwest Alentejo and Algarve coastlines [15,30]. 
Conversely, the Lusitania endemic L. nydeggeri and 
L. lanceolatum have more restricted distributions; the 
former is restricted to West and Southwest Atlantic 
sea-cliffs whereas the latter is found in the Southwest 
and South coastlines [28,30]. Limonium tetraploid taxa 
include among others, the Lusitania endemic apomict, 
L. multiflorum Erben [14,15] which exhibits both tetraploid 
and aneuploid tetraploid cytotypes (2n = 4x = 35 - [14]; 
2n = 4x = 32, 34, 35, 36 - [31]) and the aneuploid tetra- 
ploid apomict L. dodartii (Girad) O. Kuntze (2n = 4x = 35) 
which most frequently grow on maritime cliffs in the 
province of Estremadura. A third tetraploid, L. vulgare 
(2n = 4x = 36), a sexual species [19,32,33], grows in salt 
marshes [15,30]. 

Dominant genetic markers, such as those generated 
by Amplified Fragment Length Polymorphism (AFLP), 
are valuable for assessing genetic diversity within and 
between populations [34] and for inferring taxon differ- 
entiation [35], especially in species for which codomi- 
nant markers are unavailable. Some recent publications 
have added data on natural epigenetic variation in ani- 
mal and plant species by sampling cytosine methylation 
using Methylation-Sensitive Amplification Polymorphism 
(MSAP) technology [36-38], a modification to the original 
AFLP protocol that compares product profiles gener- 
ated by methylation-sensitive/insensitive isoschizomers. 
Central to the technique is the differential behavior of 
the two isoschizomer restriction enzymes (Hpall and 
Mspl) in the presence of cytosine methylation in the 
CCGG context. Hpall is inactive if one or both cytosines 
are methylated at both DNA strands, but cleaves when 
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one or both cytosines are methylated in only one strand. 
Mspl, by contrast, cleaves C5mCGG but not 5mCCGG. 
Comparison of the profiles generated by each enzyme 
from each individual allows the assessment of the methy- 
lation state of the restriction sites and so provides a rela- 
tive comparison of genetic and epigenetic variability. 
Several reports suggest that only methylation marks in 
the CG context can be transmitted between generations 
and so have potential for stable adaptive significance 
[39]. Furthermore, recent data shows that although non- 
CG methylation can be inherited, only inherited CG 
methylation is inversely correlated with gene expression 
[40] . Since only Hpall is affected by methylation of this 
kind, for simplicity, in this study we refer to profiles 
from this enzyme as epigenetic (meaning potentially 
transgenerationally stable and epigenetic) whereas Mspl 
is insensitive in this sense and so can only detect trans- 
generationally relevant genetic variation. The focus of 
this study was to therefore to use MSAP analysis as the 
primary tool in comparing the extent to which genetic 
and epigenetic diversity in natural populations of dip- 
loid and tetraploid Limonium species correlate with 
species identity and ploidy. 

Results 

Morphological differentiation between diploid and 
tetraploid species 

Herbarium specimens of the diploid species L. ovalifolium, 
L. nydeggeri and L. lanceolatum and of the tetraploid 



species L. dodartii, L. multiflorum and L. vulgar e ob- 
tained from individuals sampled in natural populations 
were used for morphometric measurements. Ten diag- 
nostic characters were selected based on an exhaustive 
review of Limonium species in Southwest Europe by Erben 
[14,15], and on previous biometric studies in the Limonium 
genus [33,41,42]. 

Only one of the ten morphological variables measured 
from representatives of the six species fitted a normal 
distribution, Maximum inner bract length (MIBL), whilst 
the other nine failed to do so, even after a logarithmic 
transformation. The remaining analyses were therefore 
performed using the original (untransformed) values 
(Additional file 1). Canonical discriminant analysis (CDA) 
of the morphological variables accounted for most of the 
variation (82.2% in the first two dimensions, comprising 
61.9 and 20.3%, respectively) and correctly assigned indi- 
viduals to species in 92.7% of cases (n=110; Figure 1; 
Table 1). However, a small number of intermediate or 
ambivalent specimens were encountered (Wilks lambda = 
0.002, x 2 = 650.259 >Xo.o5;5o = 67.505, P< 0.001). The first 
axis distinguished diploid from tetraploid species through 
the following characters: Maximum outer bract length 
(MOBL), Maximum calyx length (MCL), Maximum 
middle bract length (MMBL) and MIBL (Table 2; 
Figure 2). The second axis separated L. nydeggeri from 
the other diploid species and L. multiflorum from the 
other tetraploid species by the following characters: Max- 
imum middle bract width (MMBW), Maximum outer 
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Figure 1 Discriminant function analyses of morphometric data with predefined diploid and tetraploid Limonium species. Individuals from 
each diploid (L lanceolatum, L nydeggeri, L ovalifolium) or tetraploid (L dodartii, L multiflorum and L vulgare) species are represented by colour 
symbols. Each species centroid is represented by filled squares. Percentages of total variance explained by the functions are given in parentheses. 
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Table 1 Summary of the discriminant analysis of six predefined diploid and tetraploid Limonium species 



Species 



Predicted group membership-classification results 



Count 



Original 



L lanceolatum L nydeggeri L ovalifolium L dodartii L multiflorum L vulgare 



L lanceolatum 
L nydeggeri 
L ovalifolium 
L dodartii 
L multiflorum 
L vulgare 



0 

20 
0 
0 
0 
0 



0 
0 
0 
16 
1 
1 



0 
0 
0 
0 

2 

20 



Total 

11 

20 
21 
16 
21 
21 



% 



L lanceolatum 
L nydeggeri 
L ovalifolium 
L dodartii 
L multiflorum 
L vulgare 



90.9 

0 
14.3 

0 

0 

0 



0 

100 
0 
0 
0 
0 



9.1 

0 
85.7 
0 
0 
0 



0 
0 
0 

100 
4.8 
4.8 



0 
0 
0 
0 
85.7 
0 



0 
0 
0 
0 

9.5 
95.2 



100 
100 
100 
100 
100 
100 



Individuals from each diploid (L lanceolatum, L nydeggeri, L ovalifolium) or tetraploid (L dodartii, L multiflorum and L vulgare) species are classified. 



bract width (MOBW) and Maximum inner bract width 
(MIBW). Thus, these features were largely responsible 
for separation of species sharing the same ploidy level 

MSAP profiles analysis 

Genetic-epigenetic analyses were performed on 125 indi- 
viduals selected from natural populations within desig- 
nated SCI(s) (Figure 3; Table 3). The two MSAP primer 
combinations applied to all samples yielded 835 scorable 
fragments comprising 792 from Mspl and 778 from 

Table 2 Pooled within-groups correlations between 
discriminating variables and standardized canonical 
discriminant functions of morphological characters 



Function 





1 


2 


3 


4 


5 


MOBL 1 


0.586* 


0.414 


-0.047 


0.37 


0.311 


MCL 


0.569* 


-0.44 


-0.203 


-0.165 


-0.123 


MMBL 


0.512* 


0.298 


0.336 


-0.434 


0.032 


MIBL 


0.350* 


0.144 


-0.203 


0.193 


0.021 


MMBW 


0.405 


0.537* 


0.306 


-0.348 


0.151 


MOBW 


0.256 


0.527* 


-0.045 


0.318 


0.343 


MIBW 


-0.02 


0.442* 


0.01 


0.223 


-0139 


MNSC 


-0.334 


0.288 


-0.603 


-0.372 


0.418* 


MSL 


-0.304 


-0.08 


0.39 


0.195 


0.656* 


MNFS 


-0.005 


0.174 


0.117 


0.075 


-0.493 



Variables were ordered by absolute size of correlation within function. 
^SL, Maximum spike length; MNSC, Maximum number of spikelets per cm; 
MNFS, Maximum number of florets per spikelet; MOBL, Maximum outer bract 
length; MOBW, Maximum outer bract width; MMBL, Maximum middle bract 
length; MMBW, Maximum middle bract width; MIBW, Maximum inner bract 
width; MCL, Maximum calyx length. 

*Largest absolute correlation between each variable in respective function. 



Hpall, 92.78% and 95.36% respectively, were polymorphic 
(i.e. not present in all the analysed samples/replicates 
when restricted with one of the isoschizomers). Overall 
reproducibility between biological replicates was 83% 
and 85% for primer combinations El/Hl and E1/H3 
respectively. The methylation insensitive (genetic vari- 
ation) profiles showed only very slightly higher con- 
cordance among replicates (83.6% and 85.7%) than did 
the methylation-sensitive (epigenetic variation) profiles 
(82.4% and 84.3%). Technical reproducibility of the 
MSAP technique revealed between 92-95% band con- 
cordances (data not shown), indicating that the higher 
variability between independent DNA extractions prob- 
ably arises from variation between tissues and cell mix- 
tures used in the DNA extraction. 

Profiles from the tetraploid species (i.e. L. multiflorum, 
L. dodartii, L. vulgare) included a higher number of 
Mspl fragments (genetic profiles) per individual for both 
primer combinations (Table 4). Conversely, these species 
contained a lower number of T^^II-generated (methylation- 
sensitive epigenetic profiles) and fewer fragments per in- 
dividual than the three diploid species (i.e. L. nydeggeri, 
L. ovalifolium, L. lanceolatum), implying a higher level 
of genome-wide methylation among tetraploids. 

Genetic/epigenetic divergence of diploid and 
tetraploid species 

Principal Coordinate Euclidean Analysis (PCoA) was used 
to provide an overview of the genetic/epigenetic variability 
and structure of the studied taxa. Overall, epigenetic 
profiles {Hpall) provided imperfect but slightly better 
separation of the taxa than did the genetic profiles 
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Figure 2 Box plots of significant characters that discriminate diploid from tetraploid Limonium species. The box from diploids L nydeggerl, 
L ovolifolium, L lonceolotum, and tetraploids L dodartii, L multiflorum, L vulgore show the twenty-fifth and seventy-fifth percentile ranges and the 
median; circles and asterisks are outliers (cases with values between 1.5 and 3 box lengths from the upper or lower edge of the box). MOBL - 
Maximum outer bract length; MCL - Maximum calyx length; MMBL - Maximum middle bract length; and MIBL - Maximum inner bract length. 



(Mspl) (Figure 4A; Additional file 2). Genetic data (Mspl) 
alone also largely failed to discriminate within the diploid 
species and the tetraploid species (Figure 4B; Additional 
file 2). More remarkably, PCoA plots of the epigenetic 
Hpall profiles revealed clear separation between the 
diploid and tetraploid taxa (Figure 4C; Additional file 2). 

On average, genetic and epigenetic distances between 
the three diploid species were significantly lower than 
comparable distances between the tetraploid species 
(T-test, Genetic distance p < 0.023; epigenetic distance 
p < 0.0006 (two-tailed test); Additional file 3). Further- 
more, calculated epigenetic distances between individuals 
from different diploid species and between tetraploid spe- 
cies were higher than genetic distances between the same 
pairings (T-test, p < 0.04 (between diploid species); 



p<0.15 (between tetraploid species) (two-tailed test); 
Additional file 3). 

Analysis of genetic/epigenetic variability using Analysis 
of Molecular Variance (AMOVA) showed that 24% of 
the total observed variability can be attributed to differences 
between epigenetic {Hpall) and genetic {Mspl) sources of 
variability (Figure 4 A; Additional file 2). Independent 
analysis of results generated by each enzyme type re- 
vealed that differences between diploid and tetraploid 
samples accounted for 6-8% of the genetic variability 
and for 3-6% of the epigenetic variability. This compared 
with a much higher level of within-species variation, 
which comprised 79-83% and 74-78% of the genetic and 
epigenetic variability respectively (Figure 4B-C; Additional 
file 2). Surprisingly, while genetic differences between 
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Figure 3 Diploid and tetraploid Limonium populations sampled in continental Portugal in Sites of Community Interest. Diploid species 
(L lanceolatum, L nydeggeri and L ovalifolium) are represented by full stars and tetraploid species (L dodartii, L multiflorum and L vulgare) are 
represented by full circles. Specimens were sampled in Estremadura (PTCON0056: Baleal, Papoa, Nossa Sra dos Remedios, Vale dos Frades, Foz do 
Lizandro; PTCON0008: Cabo Raso), Alentejo (PTCON0012: Vila Nova de Milfontes, Pontal da Carrapateira), and in Algarve (PTCON0012: Cabo de 
Sagres; Praia da Luz) provinces. 



Table 3 Collection data of fourteen populations of Limonium species included in the MSAP analyses 



Ploidy level 


Species 


Site location collector* 


Geographical coordinates latitude/longitude 


N 


Diploid 


L lanceolatum 


Odemira, Vila Nova de MilFontes, ADC, APP, ASR 


37,727756/-8,770931 


10 






Peniche, llha do Baleal, ADC, APP, ASR 


39,37891 9/-9,340983 


9 






Peniche, Nossa Sr a dos Remedios, ADC, APP, ASR 


39,369906/-9,395731 


7 




L nydeggeri 


Cascais, Cabo Raso, ADC, APP, ASR 


38,71 0039/-9,485883 


7 






Aljezur, Pontal de Carrapateira, ADC, APP, ASR 


37,1950397-8,91 1 103 


10 






Vila do Bispo, Cabo de Sao Vicente, ADC ASR 


37,00261 1/-8,996564 


10 




L ovalifolium 


Vila do Bispo, Cabo de Sagres, ADC, APP, ASR 


36,994242/-8,948756 


6 






Lagos, Praia da Luz ADC, ASR 


37,087442/-8729094 


10 


Tetraploid 


L dodartii 


Odemira, Cabo Sardao, ADC, APP, ASR 


37,598344/-8,818272 


10 




L multiflorum 


Peniche, Peninsula da Papoa, ADC, APP, ASR 


39,3741 3 1/-9,377428 


6 






Lourinha, Vale dos Frades, ADC, APP, ASR 


39,276506/-9,335839 


7 






Mafra, Foz do Lizandro, ADC, APP, ASR 


38,941 53 1/-9,41 5223 


9 






Cascais, Cabo Raso, ADC, APP, ASR 


38,71 0039/-9,485883 


14 




L vulgare 


Odemira, Vila Nova Mil Fontes, ADC, APP, ASR 


37,7277567-8770931 


10 



The geographical location of each population is represented in Figure 3. Geographical coordinates of each population and sampling size (A/; approximately 10 
individuals per population) are included. 

^Abbreviations of collectors: ADC, AD Caperta: APP AP Paes; ASR, AS Rois; AC, Ana Cortinhas. Centro de Botanica Aplicada Agricultura, Instituto Superior de 
Agronomia, Lisboa, Portugal. 
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Table 4 MSAP fragment number analysis 



PlniHv IpvpI 

r iuiu y iwvwi 


Species 


E1/H1 


E1/H3 


Total 






HpaW 


Msp\ 


HpaW 


Msp\ 


Hpall Msp\ 


Diploid 


L lanceolatum 


52 


52 


47 


38 


99 90 




L nydeggeri 


50 


48 


43 


33 


93 81 




L ovalifolium 


48 


43 


42 


28 


90 71 


Tetraploid 


L dodartii 


37 


77 


24 


56 


61 133 




L multiflorum 


33 


59 


30 


35 


63 94 




L vulgare 


26 


46 


20 


39 


46 85 



Average number of MSAP fragments per diploid (L lanceolatum, L nydeggeri, 
L. ovalifolium) and tetraploid (L dodartii, L. multiflorum and L. vulgare) species 
obtained using isoschizomers enzymes HpaW (methylation sensitive) and Msp\ 
(methylation insensitive) and primer combinations E1/H1 and E1/H3. 



species accounted for 9-15%, epigenetic differences 
accounted for 19-20%, suggesting that characterisation 
of these very closely related species is best served by 
considering both genetic and epigenetic information 
rather than genetic information only (Additional file 3). 

Tetraploid species compensate lower genetic variability 
with higher epigenetic variability 

Closer consideration of the genetic and epigenetic vari- 
ability revealed that all species had more epigenetic than 
genetic variability. And, the diploid species were more 
variable (both at genetic and epigenetic level) than were 
the tetraploid ones (Table 5). However, the difference 
between genetic and epigenetic variability was higher 
among the tetraploid species despite containing lower 
levels of variability overall (Table 5). 

PCoA of the genetic profiles revealed no or only weak 
co-clustering of individuals according to population ori- 
gin, with the best example being seen among L. nydeggeri 
samples (Figure 5 A-B). No co-clustering was evident 
from the epigenetic profiles (data not shown). A lack of 
structuring according to population origin was further 
supported by AMOVA, with the main component of 
both genetic (90-97%) and epigenetic variance (86-95%) 
residing within populations. 

We next sought structuring across a geographic scale. 
Mantel test analysis revealed a correlation between genetic 
distances and geographic separation among conspecific 
populations for L. nydeggeri (R 2 = 0.784, P < 0.03). This 
correlation was significant for both primer pairs (HI/ 
El, R 2 = 0.858, P<0.02; H1/E3, R 2 = 0.616, P < 0.04) 
(Figure 6 A-B). Both primers generated a scatterplot 
showing a positive and monotonic relationship over 
all geographic distances of separation. In contrast, 
L. multiflorum populations showed no detectable rela- 
tionship between genetic distance and geographic dis- 
tance and large variance in estimates of divergence. None 
of the studied species showed a significant correlation 
between epigenetic and geographic distances, with exten- 
sive scatter between the plotted samples (data not shown). 



Finally, Mantel test analysis revealed a strong signifi- 
cant positive correlation between genetic and epigenetic 
distances among L. nydeggeri populations (Hl/El, R 2 = 
0.800, P < 0.008; H1/E3, R2 = 0.200, P < 0.07 and R 2 = 0.600, 
P < 0.005, using the information provided by both pri- 
mer combinations) (Figure 6 C). Conversely, there was 
only a weak negative correlation between genetic and 
epigenetic distances among L. multiflorum populations 
(data not shown). 

Discussion 

The effects of hybridization, polyploidy and apomixis 
have all combined to shape radiation currently seen in 
Limonium species [26,31,43]. This evolutionary model 
has been used to explain multiple series of complex ag- 
gregates of sexual diploid species and asexual polyploid 
hybrids which are perpetuated through gametophytic apo- 
mixis [44-48]. In other plant groups, hybridization and 
polyploidy have combined to generate genetic and pheno- 
typic complexity in the form of classical polyploidy pillar 
complexes. In these situations species delimitations typic- 
ally become blurred at the higher ploidy levels and inter- 
ploidy discrimination can also become difficult for some 
taxa. The many examples of this type of species complex 
include the Festuca ovina aggregate [49], Dactylis glomer- 
ata [50] and Knautia arvensis species groups [51]. The 
additional and intermittent appearance of facultative or 
obligatory apomixis, as seen in triploid and tetraploid 
Limonium species [20,21], adds another layer of com- 
plexity for species delimitation and diagnosis. In these 
instances the phenotypic range of taxa can be highly 
variable, as can their morphological distinctiveness and 
stability (e.g. [52]). The cumulative effect of these pro- 
cesses is typically manifest in the recognition of a large 
number of microspecies, with the diagnosis and genetic 
characterisation of many of the resultant taxa often being 
highly demanding (e.g. Limonium, reviewed in [26]). The 
publication of a comprehensive revision of Limonium 
species in Southwest Europe by Erben [14] in which 59 
species (13 new to science) were fully described allowed 
for several taxometric studies of geographically related 
species in this genus to be conducted. For example, vari- 
ous authors used a taxometric approach to study several 
sexual and agamospermous species of the L. binervosum 
(G. E. Sm.) C. E. Salmon complex from Western Europe 
[53], and also to examine the closely related L. vulgare 
and L. humile species in the British Isles [33]. Morpho- 
metry studies have also been used at population level. 
For instance, some authors studied the taxometric re- 
lationships between triploid or aneuploid tetraploid 
L. binervosum agamospermous colonies in the British 
Isles [42,53]. A similar strategy is applied to describe 
morphological differentiation patterns between individuals 
of a L. dufourii population from Eastern Spain [54]. 
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(See figure on previous page.) 

Figure 4 Principal Coordinate Analysis (PCoA) representing genetic and epigenetic variability in diploid and tetraploid Limonium species. 

PCoA was based on presence/absence scores of 488 polymorphic loci obtained from MSAP profiles using isoschizomers Msp\ (methylation insensitive - 
red symbols in A and B) or HpaW (methylation sensitive - blue symbols in A and C) as frequent cutters and amplified with primers (E1/H3). The first 
two coordinates were extracted and plotted against each other. Percentage of the variability shown by each coordinate is indicated between 
parentheses. Diploid species are represented by solid symbols (L lonceolotum, triangles; L. nydeggeri, rhomboids; L ovolifolium, rectangles) and 
tetraploid species are represented by empty symbols (L dodartii, triangles; L multiflorum, rectangles; L vulgare, circles). 



In the current study, CDA was applied using ten 
morphometric traits collected from representative in- 
dividuals of three diploid and three tetraploid Limonium 
species from Portugal As previous studies had suggested 
[14,33,41,42], the collective use of these characters in 
CDA was sufficient to provide clear morphological dif- 
ferentiation between species at each ploidy level This 
differentiation is based primarily on the use of seven 
morphometric variables, viz: MOBL, MCL, MMBL and 
MIBL for the first axis and MMBW, MOBW and MIBW 
on the second. Overall, these analyses not only provided 
clear separation of all diploid species but also indicated 
that L. lanceolatum and L. ovalifolium share a closer 
phenotypic affinity. Separation of the tetraploid species 
was also possible but with intermediate individuals blurring 



the boundaries between L. vulgare and L. dodartii, and 
between the latter and L. multiflorum. The rather sur- 
prising finding, however, lay in the clear phenotypic 
separation of tetraploid species from the diploids on 
the basis of the bract and calyx characteristics MCL, 
MOBL, MMBL and MIBL rather than plant size fea- 
tures more usually associated with ploidy level changes 
(e.g. see [55]). The relative importance of genetic and 
epigenetic processes in shaping the observed phenotype 
structuring among species and ploidy levels in Limonium 
has thus far remained elusive. 

Previous works on the analysis of genetic variation and 
population genetic structure in other Limonium species 
have deployed a wide range of molecular marker systems 
to infer the importance of genetic structuring in defining 



Table 5 Genetic and epigenetic variability within diploid and tetraploid Limonium species 
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f p*q) for diploid (L lanceolatum, L nydeggeri, 
(v.6.4). 
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Figure 5 Principal Coordinate Analysis (PCoA) representing genetic variability in Limonium nydeggeri populations. PCoA was based on 
presence/absence scores of 488 (primer I, El /HI) (A) and 347 (primer II, E1/H3) (B) polymorphic loci obtained from MSAP profiles using 
isoschizomers Msp\ (methylation insensitive) or HpaW (methylation sensitive) as frequent cutters. The first two coordinates were extracted and 
plotted against each other. 
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interspecies delimitation. There has, nevertheless, been 
considerable evidence that modest but significant levels 
of genetic variation does occur within species in the 
genus. For example, studies on the presumed agamo- 
spermous triploid L. dufourii have invariably revealed 
low but substantial inter- individual genetic variation in 
habitats with significant fragmentation and low popula- 
tion sizes [54,56-58]. Similarly, in diploid L. dendroides 
from Canary Islands, despite radical habitat fragmenta- 
tion and small population size, some subpopulations of 
have enough genetic variation to compensate for the in- 
fluence of drift [59]. However, in plant populations with 
low genetic variability, epigenetic variation can also act 
as an important source of potentially adaptive pheno- 
typic variability [11-13]. Nevertheless, the extent and im- 
portance of epigenetic variation in natural populations 
of sexual and agamospermous Limonium species is still 
largely unexplored [9,10,60]. 



Data generated in the current study unsurprisingly re- 
veals that tetraploids have a higher level of methylation 
than diploids. In part this may be accounted for the in- 
creased genome size, with scope for intergenome hetero- 
zysity for methylation marks adding to that expected 
between homologous loci. There are also additional in- 
direct processes that can contribute to an expected in- 
creased incidence in methylation across polyploid species. 
In newly formed polyploids, genomic duplications from 
transposons and duplication of regulatory genes, can re- 
sult in high levels of cytosine methylation on the new por- 
tions of the genome [61-63]. Several methylation-based 
processes have been implicated in the competitive silen- 
cing of duplicated gene regions (e.g. [64]), leading to 
increased methylation relative to the diploid progenitor (s). 
The same processes also lead to the expectation of dif- 
ferential methylation patterning between the diploid 
progenitor and the polyploid offspring, partly because of 



Rois et al. BMC Plant Biology 201 3, 13:205 
http://www.biomedcentral.com/1471-2229/13/205 



Page 11 of 16 



300 
250 
200 
§150 
100 
50 



A 






♦ 

y=3070.2x- 11.27 






♦ 


R 2 = 0.6156 




♦ 







0.000 



0.020 



0.040 0.060 
GDMspl 1 



0.080 



0.100 




0.000 



0.010 



0.020 



0.030 
GDMspl 2 



0.040 



0.050 



0.060 



0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 
0.00 



y=0.6432x + 0.0105 
R 2 = 0.5788 



0.000 0.020 0.040 0.060 0.080 0.100 

EpiGD (Hpoll) 

Figure 6 Correlation between pairwise genetic differentiation (GD, PhiPT) and geographical distance (GGD, in Km) (A-B) and between 
pairwise genetic and epigenetic differentiation (GD and EpiGD, respectively) (C) between Limonium nydeggeri populations. Mantle tests 
were based on MSAP data obtained using Hpoll and primer combinations primer I (El /HI) (A) and primer II (E1/H3) (B) and Hpoll and Mspl and 
both primers (C). Shown equations are the linear functions and R 2 values for each Mantel test. Analysis using 1000 permutation tests showed a 
significant correlations (A: P < 0.04; B: P < 0.02; C: P < 0.03). 



imperfect duplication of genes in the neopolyploid, as 
was reported for the Waxy gene in Spartina [64], but 
also because of various systems of asymmetric homeo- 
log silencing (e.g. [64]). Given the causal links between 
methylation and gene regulation [5], there is scope for 
such changes to influence the phenotype of the poly- 
ploids relative to their diploid relatives. Principal coord- 
inate and genetic distance analyses performed in the 
current study yielded greater separation between diploid 
and tetraploid taxa when using epigenetic information 
than when using only genetic information, suggesting 
that ploidy levels are better separated using epigenetic 
information than genetic information alone. Intriguingly, 
this pattern of variation was mirrored by the more 



pronounced morphological separation between the dip- 
loid and tetraploid taxa than that seen between species 
at either ploidy level. Should this observation apply 
more broadly to other groups, it implies that epigenetic 
profiling may provide a useful additional tool to infer 
ploidy level of preserved specimens. 

We found similar trends when making comparisons 
at the species level. Moreover, the diploid species within 
L. ovalifolium complex have imperfect but reasonable 
morphological differentiation but genetic co-variation with 
species identity was relatively modest. Interspecies separ- 
ation was more strongly enhanced when analysis was 
focused on the epigenetic variation encompassed in the 
Hpall profiles, implying that epigenetic patterning (and 
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associated gene silencing) may play a significant role in 
species separation of the group and so may have some 
utility for species diagnosis. In both cases, further re- 
search would be required to convert these anonymous 
profiles into Sequence Tagged Site epigenetic markers 
for robust diagnostic purposes. 

Analysis of genetic/epigenetic variability using AMOVA 
revealed that tetraploid species present lower levels of 
genetic variability than diploid species. This observation 
could be most plausibly explained through consideration 
of their reproductive biology. Several previous studies of 
the diploids L. ovalifolium, L. nydeggeri, and the tetra- 
ploids L. dodartii, L. multiflorum and L. vulgare have 
provided some insights into their primary reproductive 
strategies [14,31,65]. These and other works were 
based on the determination of flower dimorphisms 
linked to a sporophytic self-incompatibility system 
[14,17,18,28,33,66,67]. In these studies diploids L. nydeg- 
geri and L. ovalifolium were deemed probable sexual spe- 
cies based on their reproductive characteristics. The same 
applies to the tetraploid L. vulgare [14,65]. Conversely, tet- 
raploids L. dodartii and L. multiflorum both belonging 
to the Limonium binervosum group, were considered as 
agamospermous [31,42,53]. Hence, it seems most likely 
that the lower level of genetic diversity in the putative 
agamospermous tetraploids could be best explained by 
their apomictic reproduction mode. In other polyploid 
apomictic species, such as in Ranunculus sp., genetically 
uniform populations have been similarly observed as a 
consequence of this mode of reproduction [45]. This 
finding contrasts with the relatively high level of epigen- 
etic variability among the tetraploids, leading us to specu- 
late that in apomictic polyploid Limonium species, the 
lack of genetic variability caused by the loss of meiotic 
segregation could be partially compensated by enhanced 
epigenetic variation. Several authors have suggested simi- 
lar heritable phenotypic variation due to stable epigenetic 
variation in the absence of genetic variation (e.g., Viola 
cazorlensis; Laguncularia racemosa, and in triploid asex- 
ual dandelion lineages (reviewed in [8-10])). Viewed in this 
context, the results of the present study add Limonium to 
that list. 

The present work failed to support any relationship 
between genetic or epigenetic distance with geographic 
distance between populations of each species, except for 
a positive correlation between genetic and geographic 
distances among L. nydeggeri, consistent with regional 
equilibrium between gene flow and drift [68]. The ab- 
sence of co-correlation is not unexpected for the apomic- 
tic polyploids, and can be explained by restricted gene 
flow between populations, founder events produced by a 
limited number of individuals, absence of recombination 
and spread of single asexual clones within populations 
[45]. One plausible explanation of these results considered 



collectively is to propose that genetic information flows 
between populations but that epigenetic information is 
mainly induced locally by the environment. Alternatively, 
it might be that selection pressure on epiloci is higher 
than on genetic loci, or simply that epiloci are plastic, in 
the sense that they appear and disappear depending on 
the environmental cues. Conversely, Mantel test analysis 
of the correlation between genetic and epigenetic dis- 
tances between L. nydeggeri populations, which presents 
a regional equilibrium at genetic level, show a strong 
significant positive correlation between genetic and epi- 
genetic distances. While, genetic and epigenetic distances 
between L. multiflorum populations, which does not 
present genetic regional equilibrium, showed a weak 
negative correlation. Again, this might suggest that on 
the absence of a strong gene flow between populations, 
environmental conditions exert a higher pressure on the 
fixation of epigenetic loci that cannot be masked by gen- 
etic variability introduced by sexual reproduction. 

Conclusions 

Higher correlation was found between morphometric and 
epigenetic differentiation than between morphometric and 
genetic differentiation. We therefore, suggest that epigen- 
etic variation might be a driver of the observed phenotypic 
divergence between the studied taxa through intergenome 
silencing. We argue that the present work helps to dem- 
onstrate the importance of considering phenotypic, gen- 
etic and epigenetic variables when seeking to explain 
the dynamics of complex plant groups that feature 
hybridization, polyploidy and variable modes of repro- 
ductive biology. 

Methods 

Study species 

Natural populations of the three diploid species from 
the L. ovalifolium complex (L. ovalifolium, L. nydeggeri, 
L. lanceolatum) [28,30] and three tetraploids species 
(L. dodartii, L. multiflorum, L. vulgare) [14,15] were sur- 
veyed in the three Portuguese provinces of Estremadura 
(West), Alentejo (South- West) and Algarve (South). With 
the exception of L. lanceolatum and L. vulgare (which 
grow in salt marshes), all populations vegetated limestone 
sea-cliffs, in crevices within exposed rocks or on shallow 
soil above the rock strata and on scree slopes where 
competition with other species is very low. The loca- 
tions of all populations were recorded using Global 
Positioning System. Google Earth 6.0.2 was used for 
georeferencing and to estimate geographic distances be- 
tween populations. Geographic mapping of the popula- 
tions was performed using ArcGIS Desktop 10 (ESRI). 

Three leaves per individual were sampled from approxi- 
mately ten individuals for each population, with a total of 
125 plants included in this study. Most individuals were 
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selected from populations within designated SCI(s), 
namely in SCI Peniche/Sta Cruz (PTCON0056), SCI 
Sintra/Cascais (PTCON0008) and SCI Costa Sudoeste 
(PTCON0012) [29]. In all, five L. nydeggeri populations 
were surveyed from Ilha do Baleal (Estremadura: Peniche; 
SCI PTCON0056), Nossa Sr a dos Remedios (Estremadura: 
Peniche SCI PTCON0056), Cabo Raso (Estremadura: 
Cascais; SCI PTCON0008), Pontal da Carrapateira 
(Algarve: Aljezur; SCI PTCON0012), Cabo de Sao Vicente 
(Algarve: Vila do Bispo; SCI PTCON0012); two L ovalifo- 
lium populations from Cabo de Sagres (Algarve: Vila do 
Bispo; SCI PTCON0012), and Praia da Luz (Algarve: 
Lagos); one L. lanceolatum population from Vila Nova 
de Mil Fontes (Alentejo: Odemira; SCI PTCON0012); 
one L. dodartii population from Cabo Sardao (Alentejo: 
Odemira; SCI PTCON0012); four Limonium multi- 
florum populations from Peninsula da Papoa (Estremadura: 
Peniche; SCI PTCON0056), Vale dos Frades (Estremadura: 
Lourinha; SCI PT PTCON0056), Foz do Lizandro 
(Estremadura: Mafra; SCI PTCON0056) and Cabo Raso 
(Estremadura: Cascais; SCI PTCON0008); and one 
L. vulgare population from Vila Nova de Mil Fontes 
(Alentejo: Odemira; SCI PTCON0012) (see Figure 3; 
Table 3). Three leaves from plants at the same pheno- 
logical stage were sampled from all sites during 2010 
and kept on silica gel In this way, variation in DNA 
methylation profiles attributable to developmental or 
storage conditions differences was minimized. 

Morphometric analysis 

Morphometric analyses were performed in approximately 
twenty herbarium specimens from each species deposited 
in the herbaria Joao de Carvalho e Vasconcellos (LISI; 
ISA), Portugal. These specimens were previously collected 
in the same populations selected for MSAP analysis and 
identified on the basis of species descriptions, diagnostic 
keys, and locations already described [14], and by com- 
parison with other herbarium specimens present in 
Portuguese herbaria. The following diagnostic charac- 
ters were measured: maximum spike length (MSL), 
maximum number of spikelets per cm (MNSC), max- 
imum number of florets per spikelet (MNFS), maximum 
outer bract length (MOBL), maximum outer bract width 
(MOBW), maximum middle bract length (MMBL), max- 
imum middle bract width (MMBW), maximum inner bract 
length (MIBL), maximum inner bract width (MIBW), and 
maximum calyx length (MCL). All traits were measured in 
the lab, after removal of flower parts of each individual. 
Statistical evaluations were performed with the program 
SPSS 20 (IBM SPSS, 2010) for Windows. The morphomet- 
ric variables were tested for deviations from a normal distri- 
bution using a Kolmogorov-Smirnov test, and then these 
variables were log transformed. CDA were conducted to 
give indication of the degree to which the species were 



distinguishable from each other and to determine which 
characters contributed to this discrimination. The box- 
plots showing the medians and interquartile ranges were 
produced for each significant character for each species. 

DNA isolation 

Three replicate DNA extractions from leaves of each 
sample were performed from c. 0.05 g silica gel dried leaf 
material using the using the DNeasy 96 Plant Kit (Qiagen, 
UK) and the Mixer Mill MM 300 (Retsch, Haan, Germany) 
according to the manufacturers' instructions. DNA quality 
and quantity were verified using the nanodrop 2000 spec- 
trophotometer (ThermoScientific, Wilmington, USA). Iso- 
lated DNA was diluted in nanopure water to produce 
working stocks of approximately 10 ng uT 1 . 

Genetic/Epigenetic analyses - MSAP procedure 

We used a modification of the MSAP protocol to reveal 
global variability in CG methylation patterns between 
samples of the different specimens studied [36]. For each 
individual, 50 ng of DNA was first digested and ligated 
using 5U of EcoRI and 1U of Mspl or Hpall (New England 
Biolabs), 0.45 uM EcoRI adaptor, 4.5 uM Hpall adaptor 
and 1U of T4 DNA ligase (Sigma) in 11 uL total volume 
of IX T4 DNA ligase buffer (Sigma), 1 uL of 0.5 M 
NaCl, supplemented with 0.5 uL at 1 mg/ml of BSA for 
2 h at 37°C. The enzymes were then inactivated by 
heating to 75°C for 15 min. Following restriction and 
adaptor ligation, there followed two successive rounds 
of PCR amplification. For pre-selective amplification, 
0.3 uL of the restriction/ligation products described above 
were incubated in 12.5 uL volumes containing IX Biomix 
(Bioline, London, UK) with 0.05 uL of PreampEcoRI pri- 
mer and 0.25 uL Preamp//j^zII/Ms/?I (both primers at 
10 1M) (Additional file 4) supplemented with 0.1 uL at 

1 mg/ml of BSA. PCR conditions were 2 min at 72°C 
followed by 30 cycles of 94°C for 30 s, 56°C for 30 s and 
72°C for 2 min with a final extension step of 10 min at 
72°C. Selective PCRs were then performed using 0.3 uL 
of pre-selective PCR products and the same reagents as 
deployed for the pre-selective reactions but using FAM 
labeled selective primers (Additional file 4) Cycling 
conditions for selective PCR were as follows: 94°C for 

2 min, 13 cycles of 94°C for 30 s, 65°C (decreasing by 
0.7°C each cycle) for 30 s, and 72°C for 2 min, followed 
by 24 cycles of 94°C for 30 s, 56°C for 30 s, and 72°C 
for 2 min, ending with 72°C for 10 min. Initially, eight 
selective primer combinations (Additional file 4) were 
evaluated for their ability to detect of inter-specific vari- 
ation and to generate informative and consistent MSAP 
profiles using two replicated samples from six different 
populations (data not shown). Two primer combinations 
(El/Hl and E1/H3; Additional file 4) were chosen for 
the comparative selective amplification. 
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Fluorescently labelled MSAP products were diluted 
1:10 in nanopure sterile water and 1 uL was combined 
with 1 uL of ROX/HiDi mix (50 uL ROX plus 1 ml of 
HiDi formamide, Applied Biosystems, USA). Samples 
were heat-denatured at 95°C for 3-5 min and snap- 
cooled on ice for 2 min. Samples were fractionated on 
an ABI PRISM 3100 at 3 kV for 22 s and at 15 kV for 
45 min. 

Data analysis 

MSAP profiles were analysed using Genemapper 4.0 
software (Applera Corporation, Norwalk, Connecticut, 
USA). For analysis of the genetic/epigenetic variability 
between samples revealed using MSAP, reproducible 
product peaks were scored as present (1) or absent (0) 
to form a raw data matrix. In order to minimize the 
occurrence of fragment size homoplasy [69] only frag- 
ments with lengths between 100 and 500 bp were con- 
sidered for the analysis. All monomorphic fragments and 
any fragments present/absent in all but one individual 
were considered uninformative and removed from all data 
sets [70]. Reproducibility was estimated calculating the 
proportion of dimorphic markers between the replicates 
of the selected samples [70]. For biological error rates, we 
compared paired MSAP profiles from two leaves of one 
plant of six populations using each primer combination. 
Technical reproducibility of the MSAP technique was 
assessed through the direct comparison of profiles de- 
rived from single DNA extractions from ten representa- 
tive genotypes. 

Analysis of Genetic/Epigenetic Variance 

Genetic and epigenetic similarity between tested samples 
was determined by PCoA [71] based on the MSAP profiles 
obtained from primer combinations El/Hl and E1/H3 
using Gen Alex (v.6.4). Different components of variability 
were obtained using GenAlex (6.4) software by grouping 
the samples in two different levels (Populations (i.e. a 
group of samples) and Regions (i.e. a group of popula- 
tions). To calculate distances between natural popula- 
tions, samples from the same natural population restricted 
with each enzyme were grouped by Populations (this gener- 
ated two populations from each natural population, one 
restricted with Mspl and one with Hpall) and then at 
higher level they were grouped into two Regions (all 
samples restricted with each enzyme). To calculate dis- 
tances between species, samples from the same taxa 
restricted with each enzyme were considered and grouped 
as one Population (this generated two populations from 
each original species, one restricted with Mspl and one 
with Hpall) and then at higher level they were grouped 
into two Regions (all samples restricted with each enzyme). 
We then used AMOVA [72] to evaluate the structure 
and degree of epigenetic diversity among and between 



populations, and between species. Pairwise PhiPT [72] 
comparisons (an analogue of the Fst fixation index, that 
measures differential connectivity/genetic diversity among 
populations) between samples restricted with Mspl or 
Hpall was used to infer their overall level of genetic or 
epigenetic divergence respectively. AMOVA was subse- 
quently calculated using GenAlex (v.6.4) to test the sig- 
nificance of PhiPT between populations and species [73] 
with the probability of non-differentiation (PhiPT = 0) be- 
ing estimated over 9,999 permutations. We then calcu- 
lated genetic diversity estimates (expected heterozygosity, 
Hj) and the actual genetic diversity for each of the groups 
above, by using Shannons index (Shi) [74]. Finally, pair- 
wise genetic or epigenetic and geographical distance 
(kilometres) matrices were analysed between populations 
within species and among species by means of a Mantel 
test. The level of significance was assigned after 1000 per- 
mutation tests, as implemented in Genalex 6 [75]. 

Additional files 



Additional file 1: Mean values of morphometric characters in diploid 
and tetraploid Limonium species. Diploid L. lonceolatum, L nydeggeri, 
L. ovolifolium and tetraploid L. dodartii, L. multiflorum and L. vulgare species 
are considered. 

Additional file 2: Principal Coordinate Analysis (PCoA) representing 
genetic and epigenetic variability in diploid and tetraploid 
Limonium species. PCoA was based on presence/absence scores of 347 
polymorphic loci obtained from MSAP profiles using isoschizomers Mspl 
(methylation insensitive - red symbols in A and B) or HpaW (methylation 
sensitive - blue symbols in A and C) as frequent cutters and amplified 
with primers (E1/H1). The first two coordinates were extracted and plotted 
against each other. Percentage of the variability shown by each coordinate 
is indicated between parentheses. Diploid species are represented by solid 
symbols (L lanceolatum, triangles; L. nydeggeri, rhomboids; L. ovalifolium, 
rectangles) and tetraploid species are represented by empty symbols 
(/_. dodartii, triangles; L. multiflorum, rectangles; L. vulgare, circles). 

Additional file 3: Estimated genetic (GD) and epigenetic (EpiDG) 
distances between and within diploid and tetraploid Limonium 
species. Diploid (/_. lanceolatum, L. nydeggeri, L. ovalifolium) and tetraploid 
(L dodartii, L. multiflorum and L. vulgare) species are considered. 

Additional file 4: Oligonucleotides used for MSAP analysis. Selective 
nucleotides are indicated as + XYZ in the primer code column. Enzyme 
column indicates the restriction enzyme site associated with each primer. 
* FAM labeled selective primers. 
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